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Abstract In this paper, transverse momentum spectra 
of 7r + , K + and p measured at fix event- multiplicities 
and yfs — 0.9, 2.76 and 7 TeV collision energies by 
the CMS Collabotaion are shown to fit the Tsallis- 
distribution. It is found that the power of the distri- 
bution shows a double-logarithmic dependence on the 
event-multiplicity N, while the T parameter depends 
linearly on N. A similar double-logarithmic dependence 
of the q parameter of ir° spectra on the collision energy 
y/s is found too. 

It is also shown that event- by- event fluctuations of 
the multiplicity N and the total Et energy going into 
the transverse region can be the reason for the emer- 
gence of the Tsallis distribution in high-energy proton- 
proton collisions. 

Keywords hadron spectra • Tsallis-statistics • super- 
statistics • KNO scaling • multiplicity fluctuations • 
transverse energy fluctuations 



1 Introduction 

Cut-power law hadron distributions (sometimes refered 
to as Tsallis distribution) are observed in various high- 
energy collisions: in electron-positron e + e~ [T] - [7], 
proton- (anti-)proton pp [5] - [H], elastic pp Q25] and 
nucleus-nucleus A A [20] - [33] collisions. While the de- 
pendence of the parameters of this distribution on the 
collision energy yfi [ip||5PP2p3P5PDpip5PBP2] , 
centrality [25,26,28,33 , number of colliding nucleons 
Neon |2g], hadron species $SKMBHEMMMnMflBI\ 
and momentum transfer [19) is thoroughly discussed in 
the literature, the dependence of hadron spectra on the 
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event-multiplcity measured in [33] has not been ana- 
lysed yet. In Sec. [3] this task is fulfilled. 

One explanation for the emergence of cut-power laws 
may be found in thermal hadronisation models. These 
models are based on the conjecture that in high-energy 
collisions, small thermal droplets of matter (often ref- 
ered to as "fireballs" or "clusters") are created and 
these droplets fragment into hadrons. The calculations 
are carried through either in the canonical fT3"lll41l3"Tll3"8l 
142 1 143 ] or microcanonical [4lll ^ [lT l [34 l [35 l [36 l l4T] frame- 
work, and describe measured data on hadron specra, 
total hadron multiplicities as well as multiplicity distri- 
butions. The latter can be approximated by either the 
negative binomial or by Euler's gamma-distribution. 

Since each formed fireball is different, they may carry 
different four-momenta, fragment into different number 
of hadrons and - if one works in the canonical frame- 
work - they may have different temperature. It has been 
shown that, if a single cluster or fireball is described by 
the Boltzmann-Gibbs distribution, however, the tem- 
perature [smiiimminignigg] or cluster volume esi 

[54] or the multiplicity of hadrons stemming from a clus- 
ter [SlfTB] fluctuates from cluster to cluster, the average 
hadron spectrum may become a cut-power law (or Tsal- 
lis) distribution. The case of volume [34] and multiplic- 
ity [5ll6lH8] fluctuations has also been generalised to the 
microcanonical ensemble. The latter has been applied 
to the fragmentation of jets produced in e + e~ [5][5] and 
pp collisions [T5] . 

It is worth noting that if we approximate very nar- 
row jets with one-dimensional bunches of massless had- 
rons, the conservation of four-momentum reduces to 
the conservation of energy, and the clusters (jets) be- 
come massless too (for a discussion, see [H]). Thus, jet- 
fragmentation models [5,6,18 may be considered as an 
approximation of microcanonical statistical hadronisa- 
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tion models Pl[THlPTl|Ml[^l[Ml|iT] for very high hadron 
energies and jet-like, very narrow, one-dimensional clus- 
ters. The latter models are sensitive to the species of 
hadrons and multiplicity distributions are obtained from 
them as a consequence. In the previous models, the mul- 
tiplicity distribution is put in by hand but, because of 
their simplicity, analitic calculations are possible. 

In this paper, I will focus on the effect of clus- 
ter energy fluctuations. Hadronic transverse momen- 
tum spectra in pp collisions are mixtures of particle 
yields comming from many clusters of different energy 
E c , four- velocity u M , and yielding different N number 
of hadrons. For the spectra published in [44] is soft, 
I will work in the canonical ensemble, supposing that 
each cluster produces hadrons according to an isotrope 
Boltzmann-Gibbs distribution in D = 3 dimensions in 
the frame co-moving with the cluster: 



Gibbs distribution in the frame co-moving with the 
cluster: 



dN 



d D p 



1 cluster 



> DN 



(2) 



with E c and N being the energy and multiplicity of the 
cluster and A = [DN/E C ] U /[ KD r(D)} comes from the 
normalisation condition 



d D pf N ,E^) = 1, 



(3) 



where kd — J dfijj is the angular part of the D di- 
mensional momentumspace integral. The masses of the 
hadrons have been neglected, taking e= |p|. Let us also 
conjecture cluster energy fluctuations of the form 



dN 
d?p~ 



exp{-/3(w M p p - to)} , 



(1) 9n(E c ) 



1 
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1 cluster 



with f3 being the inverse temperature and p M being the 
four-momentum of the hadron. 

In Sec. [5J I show that an appropriate choice for 
fluctuations of the cluster energy result in cut-power 
law (or Tsallis) shaped hadron spectrum. As a first ap- 
proximation, I will neglect the cluster velocities choos- 
ing m m = (1,0,0,0) and do calculations for massless 
(m = 0) particles. This way, (3 = DN/E C holds. Fur- 
theremore, I will assume that only one cluster is formed 
in one pp event. 

In Sec. |31 I show that the Tsallis distribution ob- 
tained in Sec. [5] fits spectra of identified tt + , K + and p 
at fixed event multiplicity at LHC energies. 

In Sec. 2J I show that multiplicity fluctuations of the 
Euler-gamma type superimposed over cluster energy 
fluctuations also result in an approximately Tsallis- 
shaped transverse n + spectrum and fit measurements 
well. In addition, I also show that a double-logarithmic 
dependence of the q parameter of 7r° spectra on the col- 
lision energy y/s holds in the y/s € [0.2, 7] TeV range. 

In Sec. [SJ I make a prediction on the joint distri- 
bution of the event multiplicity N and total transverse 
energy Et in an event, supposing that Et oc E c . 



2 Power-law Spectrum from Cluster Energy 
Fluctuations 

According to Eq. (p}, let us conjecture that hadrons 
stemming from a cluster have an isotropic Boltzmann- 
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This way, the hadron distribution in a cluster of 
multiplicity N , averaged over fluctuations of the cluster 
energy, becomes a cut-power law (or Tsallis) distribu- 
tion: 



dN 



d D p 



1 cluster / E, 



= f N (e)= I dE c g N (E c )f NtEc (e) 

-1/(9-1) 



(5) 



with 



q = i + 



T = 



I 



a + D + 1 
aEo 



DN(a + D + l) ' ^ 

The parameters a and Eq (and thus q and T) may 
depend on the hadron multiplicity N in the cluster. 



3 Hadron Spectra at Fixed Multiplicity 

In this section, I show that Eq. ([5J (in D = 3 di- 
mensions) fits transverse momentum spectra of various 
identified hadrons if the kinetic energy, e — > m<r — m, is 
used as scaling variable: 



dN 



dp T dy 



Apt rriT 



[l+2^(m T -m)] 



1/(9-1) ' 



(7) 
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The analysed data are: transverse momentum spectra 
ol7r + (Fig.[TJ K + (Fig. [3]) andp (Fig. [5]) stemming from 
pp collision events of fixed multiplicities (indicated as 
Ntracks in the figures) at yfs — 0.9, 2.76 and 7 TeV 
collision energies and from the rapidity range \y\ < 1. 
Figs. HHi] and [7] show that a 

q = 1 + /j,lnla(N/N q ) , 

T = T {1 + N/N T ) (8) 

dependence of the q and T parameters on the event- 
multiplicity N is consistent with measurements. The 
parameters fx, N q , T Q and Nt take different values for 
each particle species, but do not seem to change signifi- 
cantly (within errors) as the collision energy yfs varies. 
It is interesting that while the T parameters of K + and 
p grow with TV, the T of 7r + is approximatly indepen- 
dent of the event-multiplicity. Since the dominant part 
of the produced hadrons is pions, we may use Eq. {8j| 
to estimate the dependence of the q and T parameters 
of pions on the pion- multiplicity N n , by substituting 
N — > N-x . However, it is perhaps not true for kaons and 
protons. 



pions. The analysis of transverse momentum spectra of 
pions in Sec. [3] show that the q and T parameters of 
the spectrum Eq. ([7]) depend weakly on the event mul- 
tiplicity (T ss 70 - 80 MeV almost independently of N 
and n « 0.13 - 0.14 in Eq. ©). Thus, the multiplicity 
and cluster energy averaged spectrum Eq. (|10[) , evalu- 
ated with parameters obtained from the analysis of pion 
spectra in Sec. [3J is close to the Tsallis distribution. 

Fig. [0] shows measured and calculated multiplicity- 
averaged 7r + spectra. In the evaluation of Eq. (ITU1) , the 
discrete sum was approximated by an integral with a 
lower bound 2V_ = e N q in order to ensure that q > 1. 
For the a parameter of the multiplicity distribution 
Eq. ([9j, I used values obtained from fits to experimental 
data shown in Fig. [5] For the mean multiplicity parame- 
ter No, values higher then the ones obtained from fits to 
measured multiplicity distributions turned out to give 
best agreement with measured spectra (see the caption 

of Fig. p. 

It is also worth to note that the q parameter of the 
transverse momentum spectra of neutral pions (which 
is also an E c and N averaged observable) also show 
double-logarithmic dependence on the collision energy 
for y/s € [0.2, 7] TeV 



4 The Effect of Multiplicity Fluctuations 



q{s) = l + qi lnln(Vs/Q ) 



(11) 



It has been known for a time that multiplicity distribu- 
tions of charged hadrons in high-energy collisions may 
be approximated by the negative binomial (NBD) or 
Euler' Gamma distribution [4l[51[T^[T51|3"4l|3'51[3"8l[3"9l|i0l . 
Fig. [5] shows fits of Euler's Gamma-distribution 



p(N) 



r(a) N \ N J 



aNY 



-aN/No 



(9) 



to multiplicity distributions measured by the ALICE 
Collaboration (15114?] in the \rj\ < 1 rapidity interval. 
(Eq. © is the n > 1 and k/{n) » 1 limit of the NBD 
used in [5011551 to analyse data in gHES].) 

The hadron distribution in a cluster, averaged over 
the fluctuations of both the cluster energy E c and the 
cluster multiplicity N, is 



dN 



d D p 



1 cluster/ E C ,N 



f(e)=X2 P (N)f N (e) 



(10) 



N 



For pions, the m — choice used in Sec. [5] may be 
an acceptible approximation in the measured pr range 
in |44j . Furtheremore, the multiplicity distribution of 
pions may take a similar form as the experimentally ac- 
cessible multiplicity distributions of charged particles, 
since most of the hadrons produced in pp collisions are 



In the meanwhile, the T parameter is not affected signi- 
ficantly (within errors) by the grows of s. Fig. [TO] shows 
fit results. Similar functional forms for q(s) have also 
been proposed in pllT4ll2Tll22l . 

From the observation that in pp collisions at t/s = 
900 GeV the q parameter of various hadron species coin- 
cide within errors [9 , 10 , 1 1 , we may infere that Eq. ([IT]) 
might hold for other hadons too. 



5 Joint Distribution of the Multiplicity and the 
Total Transverse Energy in an Event 

In the model presented above, the joint distribution of 
cluster energy E c and cluster multiplicity N is 



p(N,E c )=p(N)g N (E c ) 



(12) 



where p(N) is the E c - averaged multiplicity distribu- 
tion 



p(N) 



dE c p{N,E c ) 



(13) 



shown in Fig. [5] and gN{E c ) is the normalised distri- 
bution of E c at fix multiplicity. From Eqs. (JBJ) and ©, 

the parameters of gN(E c ) depend on N as 
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Hln\n{N/N q ) 

DNTp (1 + N/N T ) 
° l-(D + l)nlnln(N/N q ) ' 1 ' 

Since the hadron distribution in a cluster is assumed 
to be isotropic, the Et energy going into the transverse 
region (\y \ < 1 in |44j ) in a sigle pp event is proportional 
to the energy of the cluster formed in the event: Et ~ 
E c . Thus, the distributions of Et and E c have the same 
form. Fig. [TT1 shows the joint distribution p(N, Et) with 
its projections, using parameters obtained from data on 
7r + spectra measured at y/s = 7 TeV and in the rapidity 
range \y\ < 1. 

6 Conclusions 

Cut-power law hadron distributions (sometimes refered 
to as Tsallis distribution) are observed in various high- 
energy collisions (from electron-positron to nucleus- 
nucleus reactions [1] - [33] ) . One possible explanation 
for this phenomena can be found in thermal hadro- 
nisation models. According to these models, in high- 
energy collisions, small thermal "fireballs" or "clusters" 
are formed and fragment into hadrons. Hadrons stem- 
ming from a single cluster inherit the thermal (canon- 
ical [37U1DI1113] or microcanonical HHnETHMHSillMl 
|4"T] ) distribution of the fireball. However, since the four- 
momentum (or four- velocity) and cluster mass [H HCTTTl 
[34p5l[36ll4T]. or the temperature [rj[2l[20l[2T][22l[23l[29l 
[33] or cluster volume [231134] or the hadron multiplic- 
ity [51IT8] fluctuates from cluster to cluster, the average 
hadron spectrum may become a cut-power law (or Tsal- 
lis) distribution. 

The measurement of transverse spectra of a few 
types of hadrons stemming from proton-proton colli- 
sions for some fixed event multiplicities |44j made it 
possible to get rid of the effect of event- by- event mul- 
tiplicity fluctuations on hadron spectra. In Sec. |3l it is 
shown that transverse spectra of ir + , K + and p take 
a cut-power law shape (Eq. ([7])) even in fixed multi- 
plicity proton-proton events. It is also found that the 
q parameter of the spectra follows a double-logarithmic 
dependence on the event-multiplicity N, while the T 
parameter is independent of N (within errors) for pi- 
ons, but grows linearly with TV for kaons and protons 
(see Eq. © and Figs. 0). 

In Sec.[2j it is shown that cluster energy fluctuations 
of the form of Eq. ((4]) result in cut-power law shaped 
hadron spectrum at fix cluster multiplicity, if hadrons 
in a cluster are distributed according to the Boltzmann- 
Gibbs distribution. Averaging this spectrum over multi- 



plicity fluctuations Eq. (O , using the weak multiplicity 
dependence of the q and T parameters of pion spectra 
found in Sec. [3l also results in an approximately cut- 
power law spectrum that describes transverse tt + spec- 
trum (see Sec. [4] and Fig. [9]). This is consistent with 
the observation that multiplicity averaged transverse 
spectra (usually simply called 'transverse spectra') of 
hadrons stemming from pp collisions fit the Tsallis dis- 
tribution. In addition, I have also shown that the q 
parameter of the transverse spectrum of tt° shows a 
double-logarithmic dependence on the collision energy 
^ in the range e [0.2, 7] TeV (see Fig. [10]). Sim- 
ilar functional forms for q(s) have also been proposed 
in [T3l[T4l[2Tl[22]. 

In order to be able to perform analytic calculations, 
I conjectured that only one cluster is formed in a sin- 
gle proton-proton event; hadrons stemming from the 
cluster have an isotropic Boltzmann-Gibbs distribution 
in the frame co-moving with the cluster; finally, I ne- 
glected particle masses and cluster velocities. In this 
"first approximation" , the E c energy of the cluster for- 
med in a single proton-proton event is proportional 
to the Et energy that reaches the transverse region 
(\y\ < 1 in [M]), E c cx Et- Thus, the joint distribution 
p(N, Et) of the event multiplicity and the transverse 
energy in an event can be predicted (see Sec. [5] and 
Fig. EE]). 

Though, it is to be emphasized that there are rough 
approximations in the above presented calculations (the 
case of multiple cluster production in a proton-proton 
event as well as non-zero cluster velocities and finite 
hadron masses are to be taken into account in future 
works), it is clearly pointed out that event-by-event 
fluctuation of the hadron multiplicity as well as that 
of the energy reaching the transverse range may result 
in cut-power law (or Tsallis) shaped hadron spectra. 
Even if hadrons had Boltzmann-Gibbs distribution in 
a single event or cluster. 
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Fig. 1 Transverse momentum spectra of 7r+-s stemming Fig. 2 Dependence of the q parameter of Eq. ([7]) on the 

from pp collisions of fix multiplicities, measured at \fs = 0.9 event-multiplicity obtained from fits to 7r + spectra shown in 

TeV (top), 2.76 TeV (middle) and 7 TeV (bottom) collision Fig. CD Top: ^fs = 0.9 TeV, middle: = 2.76 TeV bottom: 

energies. Rapidity range: \y\ < 1. Data of graphs were pub- ^/s = 7 TeV. Curves are fits of Eq. JSJ. 
lished in [2]. Curves are fits of Eq. QJ. 
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Fig. 3 Transverse momentum spectra of K + -s stemming 
from pp collisions of fix multiplicities, measured at y/s = 0.9 
TeV (top), 2.76 TeV (middle) and 7 TeV (bottom) collision 
energies. Rapidity range: \y\ < 1. Data of graphs were pub- 
lished in [2]. Curves are fits of Eq. J7J. 
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Fig. 4 Dependence of the q parameter of Eq. (0 on the 
event-multiplicity obtained from fits to K + spectra shown in 
Fig. [3] Top: V« = 0.9 TeV, middle: = 2.76 TeV bottom: 
^/i = 7 TeV. Curves are fits of Eq. ®. 
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Fig. 5 Transverse momentum spectra of p-s stemming from 
pp collisions of fix multiplicities, measured at y/s = 0.9 TeV 
(top), 2.76 TeV (middle) and 7 TeV (bottom) collision ener- 
gies. Rapidity range: \y\ < 1. Data of graphs were published 
in [ii]. Curves are fits of Eq. @. 
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Fig. 6 Dependence of the q parameter of Eq. (J7]) on the 
event-multiplicity obtained from fits to p spectra shown in 
Fig. [5] Top: y/s = 0.9 TeV, middle: y/s = 2.76 TeV bottom: 
y/s = 7 TeV. Curves are fits of Eq. ®. 
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Fig. 7 Dependence of the T parameter of Eq. on the Fig. 8 Top: fits of Eq. to multiplicity distributions of 

event- multiplicity obtained from fits to 7r+ , K + and p spectra charged hadrons stemming from pp collisions at y/s = 0.9, 
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Fig. 9 Transverse momentum spectra of 7r+-s stemming 
from pp collisions and averaged over multiplicity fluctua- 
tions. Measured data are published in |44 |. Dashed lines show 
Eq. HOI) with different parametrisations of the multiplicity 
distribution Eq. J9j- For collision energies y/s = 0.9, 2.76, 
and 7 TeV, a = 1.38, 1.4, 1.12 and N = 15.3, 24, 32 were 
used respectively (see text below Eq. fllOl) ). 



4. C. Bignamini, F. Becattini, F. Piccinini, A Monte-Carlo 
generator for statistical hadronization in high energy e+e- 
collisions, Eur. Phys. J. C, 72, 2176, (2012) 

5. K. Urmossy, G. G. Barnafoldi, T. S. Biro, Generalised 
Tsallis Statistics in Electron-Positron Collisions, Phys. 
Lett. B, 701, 111-116, (2011), larXiv:1101.3023l 

6. K. Urmossy, T. S. Biro, G. G. Barnafoldi, Generalised 
Microcanonical Statistics and Fragmentation in Electron- 
Positron Collisions, Acta Physica Polonica B, 5:(2), pp. 
363-368, (2012) 

7. T.S. Biro, K. Urmossy, P. Van, G.G. Barnafoldi, Z. 
Schram, Non-extensive statistical model for strange and 
non-strange hadron spectra at RHIC and LHC energies, 
Acta Physica Polonica B, 43:(4), pp. 811-820, (2012) 

8. Cheuk-Yin Wong, G. Wilk, Tsallis Fits to p t Spectra for 
pp Collisions at LHC, Acta Physica Polonica B, 42, 2047, 
(2012), arXiv:1210.3661 

9. J. Cleymans, The Tsallis Distribution at the LHC, Pro- 
ceedings of the 1 st International Conference on New Fron- 
tiers in Physics ICFP, (2012). laxXiv: 1210.74641 

10. J. Cleymans, D. Worku, The Tsallis distribution in pro- 
tonproton collisions at y/s = 0.9 TeV at the LHC, J. Phys. 
G: Nucl. Part. Phys., 39, 025006, (2012) 

11. J. Cleymans, D. Worku, Relativistic Thermody- 
namics: Transverse Momentum Distributions in High- 
Energy P hysics, Eur. Phys. J. A, 48, 11, 160, (2012) 
larX iv: 1203.43431 

12. G. G. Barnafoldi et al., Tsallis-Pareto like distributions in 
hadron-hadron collisions, J. Phys. Conf. Ser., 270, 012008, 
(2011) 

13. T. Wibig, I. Kurp, Large Transverse Momenta in Sta- 
tistical Models of High Energy Interactions, JHEP, 0312, 
039, (2003); 

14. T. Wibig, The non-extensivity parameter of a thermody- 
namical model of hadronic interactions at LHC energies, J. 
Phys. G: Nucl. Part. Phys., 37, 115009, (2010) 



■a 



10' 



- % 


7C° 


pp @ 7000 GeV 
q = 1,129 +0.001 


r * 




T = 0.099 ± 0.003 


ft 




pp @ 900 GeV 






q = 1 .1 16 + 0.006 


-A 




T = 0.096+0.010 






pp @ 200 GeV 






q = 1.092 ±0.003 






T = 0.104 + 0.008 






pp @ 62 GeV 






q = 1.070+0.002 






T = 0.1 10 ±0.005 



20 



PT 



1.15 

cr 



q = 1 +q i lnln(Vs/Q ) 
q i = 0.063 ± 0.001 
Q„ = 2.793 + 0.406 GeV 



1.1" 



1.05- 




10 4 

Vs [GeV] 



0.2 



0.15- 



0.1 



0.05- 



T = O.I 02 + 0.002 GeV 



-i i ■ 



10* 



10 J 



Vs [GeV] 



Fig. 10 Transverse momentum spectra of 7r°-s stemming 
from pp collisions at various collision energies and fits of 
Eq. ((7J (top), y^s dependence of the q parameter and fit of 
Eq. ( jlip (middle), y/s dependence of the T parameter (bot- 
tom) . Data of graphs were published in | 45II46||4T] . 



10 



II 



10"^ 



10' 



10" : 



- 




m m m m i 


tj = Z UbV 

E T = 4 GeV 
E = 6 GeV 




/ \ 

/ V 

' • \ 
* ,* 1 


..-> 

"Ay 


E T = 8 GeV 




( - + t 
1 .' " 

f .' ' 
i; ' 


\ N '-^ 
V v 

\ \ 
\ \ 
\ % 
\ \. 






;// 


\ \ 
\ * 


S 




j ■* t 

\\\ 

w 




V *•. *" 
\ % S 
V 

V '* 

\ 

\ 



20 



40 



N 




N = 8 
N = 16 
N = 24 
N = 32 



20 



40 



E T [GeV] 



LU 



0.08 



0.06 




-0.04 



0.02 



Fig. 11 Plots of Eq. JT5J the joint distribution p(N,E T ) of 
the event-multiplicity N and total transverse energy Et (bot- 
tom) and its projections (top and middle). The parameters 
used were obtained from fits to the multiplicity distribution 
of charged hadrons and transverse 7r+ spectra measured in pp 
collisions at y/s = 7 TeV collision energy. Used parameters 
are: a = 1.11, N = 12, /j, = 0.14, N q = 1.2, T = 70 MeV, 
JV T = 5000. 



15. G. G. Barnafoldi, G. Kalmar, K. Urmossy, T. S. Biro, 
Tsallis-Pareto like distributions in hadron-hadron colli- 
sions, Gribov 80 - Memorial Volume, (ISBN: 978-981-4350- 
18-1), p. 357., World Scientific, Singapore, (2011) 

16. F. Becattini, G. Passaleva, Statistical hadronisation 
model and transverse momentum spectra of hadrons in high 
energy collisions, Eur. Phys. J. C, 23, 551-583, (2002) 

17. F. M. Liu, K. Werner, J. Aichelin, M. Bleicher, H. 
Stoecker, A Micro- Canonical Description of Hadron Pro- 
duction in Proton-Proton Collisions, J. Phys. G, 30, 589- 
594, (2004) 

18. K. Urmossy, G. G. Barnafoldi, T. S. Biro, Microcanonical 
Jet-fragmentation in proton-proton collisions at LHC En- 
ergy, Phys. Lett. B, 718, 125-129, (2012), larXiv: 1204. 15081 

19. D. A. Fagundes, M. J. Menon, P. V. R. G. Silva, Prelim- 
inary Results on the Empirical Applicability of the Tsallis 
Distribution in Elastic Hadron Scattering, arXiv:1204.5646 
(2012) 

20. M. Rybczynski, Z. Wlodarczyk, G. Wilk, On the possi- 
bility of q-scaling in high energy production processes, J. 
Phys. G: Nucl. Part. Phys., 39, 095004, (2012) 

21. G. Wilk, Z. Wlodarczyk, Consequences of temperature 
fluctuations in observables measured in high energy colli- 
sions, larXiv: 1203.4452 1 (2012) 

22. G. Wilk, Z. Wodarczyk, The imprints of superstatistics in 
multiparticle production processes, Central European Jour- 
nal of Physics, 10, Issue 3, pp. 568-575, (2011) 

23. G. Wilk, Z. Wlodarczyk, Equivalence of volume and tem- 
perature fluctuations in power-law ensembles, J. Phys. G, 
38, 065101, (2011) 

24. Z. Tang, L. Yi, L. Ruan, et al., Statistical Origin 
of Constituent-Quark Scaling in the QGP hadronization, 
larXiv:1101.19121 (2011) 

25. M. Shao, L. Yi, Z. Tang, et al., Examine the species and 
beam-energy dependence of particle spectra using Tsallis 
Statistics, J. Phys. G: 37, 085104, f2010). larXiv:0912.0993l 

26. D. D. Chinellato, J. Takahashi, I. Bediaga, A Non- 
extensive Equilibrium analysis of 7r+ pt Spectra at RHIC, 
J. Phys. G: 37, 094042, (2010), arXiv:1001.3136 

27. K. Urmossy, T. S. Biro, Cooper-Frye Formula and Non- 
extensive Coalescence at RHIC Energy, Phys. Lett. B: 689, 
14-17, (2010), larXiv:0911.14li~l 

28. B. De, G. Sau, S. K. Biswas, S. Bhattacharyya, P. Gup- 
taroy, Analyzing Non-Extensivity of -spectra in Relativistic 
Heavy Ion Collisions at ^/sjvjv = 200 GeV, J. Mod. Phys. 
A: 25, 1239-1251, (2010) 

29. G. Wilk, Z. Wlodarczyk, Power laws in elementary and 
heavy-ion collisions - A story of fluctuations and nonexten- 
sivity?, Eur. Phys. J. A, 40, 299-312, (2009) 

30. Z. Tang, Y. Xu, L. Ruan, G. Buren, F. Wang, Z. Xu, 
Spectra and radial flow in relativistic heavy-ion collisions 
with Tsallis statistics in a blast-wave description, Phys. 
Rev. C: 79, 051901 (R), (2009) 

31. T. S. Biro (ed.) et al., Eur. Phys. J. A: 40 255-344, (2009) 

32. J. Cleymans, G. Hamar, P. Levai, S. Wheaton, Near- 
thermal equilibrium with Tsallis distributions in heavy ion 
collisions, J. Phys. G: 36, 064 018, (2009) 

33. M. Biyajima, T. Mizoguchi, N. Nakajima, N. Suzuki, G. 
Wilk, Modified Hagedorn formula including temperature 
fluctuation - Estimation of temperatures at RHIC exper- 
iments, Eur. Phys. J. C: 48, 597-603, (2006) 

34. V. V. Begun, M. Gazdzicki, M. I. Gorenstein, Power Law 
in Micro-Canonical Ensemble with Scaling Volume Fluctu- 
ations, Phys. Rev. C, 78, 024904, f2008). lajXiv: 0804.0075 

35. F. Becattini, L. Ferroni, Statistical hadronization and 
hadronic microcanonical ensemble I, Eur. Phys. J. C, 35, 
243-258, (2004) 



11 



36. F. Becattini, L. Ferroni, Statistical hadronization and 
hadronic microcanonical ensemble II, Eur. Phys. J. C, 38, 
225-246, (2004) 

37. R. Hagedorn, Suppl. Nuovo Cim., 3, 147, (1965); M. 
Szczekowski, G. Wilk, Phys. Lett. B, 374, 225, (1996); K. 
Redlich, A. Andronic, F. Beutler, P. Braun-Munzinger, J. 
Stachel, J. Phys. G, 36, 064021, (2009) 

38. F. Becattini, A. Giovannini, S. Lupia, Multiplicity distri- 
butions in a thermodynamical model of hadron production 
in e+e" collisions, Z. Phys. C, 72, 491, (1996) 

39. P. Ghosh, Negative binomial multiplicity distribution in 
proton-proton collisions in limited pseudorapidity intervals 
at LHC up to t/s = 7 TeV and the clan model, Phys. Rev., 
D85, 054017, (2012) 

40. A. Dumitru, Y. Nara, KNO scaling of fluctuations in 
pp and pA, and eccentricities in heavy-ion collisions, Phys. 
Rev. C, 85, 034907, (2012) 

41. F. M. Liu, K. Werner, Micro-canonical pentaquark pro- 
duction in ee annihilations, Phys. Rev. D, 74, 034024, 
(2006) 

42. F. Becattini, U. W. Heinz, Thermal hadron production 
in p p and p anti-p collisions, Z. Phys. C, 76, 269, (1997) 

43. F. Becattini, J. Cleymans, A. Keranen, E. Suhonen, K. 
Redlich, Features of particle multiplicities and strangeness 
production in central heavy-ion collisions between 1.7A- 
GeV/c and 158A-GeV/c, Phys. Rev. C, 64, 024901, (2001) 

44. The CMS Collaboration, Study of the inclusive produc- 
tion of charged pions, kaons, and protons in pp collisions 
at sqrt(s) = 0.9, 2. 76, and 7 TeV, S ubmitted to the Eur. 
Phys. J. C (2012). fariTiv: 1207.47241 

45. Abelev et.al. [The ALICE Collaboration], Neutral pion 
and r) meson production in proton-proton collisions at \fs = 
0.9 TeV and y /s = 7 TeV, Phys. Lett. B, 717, 162, (2012), 
larXiv: 1205. 57241 

46. The PHENIX Collaboration, Inclusive cross section and 
double helicity asymmetry for piO production in p+p colli- 
sions at sqrt(s)= 62.4GeV, Phys. Rev. D, 79, 012003 (2009), 
larXiv:0810.0701l 

47. The PHENIX Collaboration, Midrapidity Neutral Pion 
Production in Proton-Proton Collisions at sqrt(s) = 200 
GeV, Phys. Rev. Lett., 91, 241803 (2003), [hep-ex/03 04038 

48. The ALICE Collaboration, Charged-particle multiplicity 
measurement in proton-proton collisions at ^Js — 7 TeV 
with ALICE at L HC, Eur. Phys. J. C, 68, 345-354, (2010), 
larXiv: 1004.35141 

49. The ALICE Collaboration, Charged-particle multiplicity 
measurement in proton-proton collisions at ^fs = 0.9 and 
2.36 TeV with ALICE at LHC, Eur. Phys. J. C, 68, 89-108, 
(2010) 



